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Abstract 

We present an unexpected application of tropical convexity to the determi- 
nation of invariants for linear systems of differential equations. We show that 
the classical Gerard-Levelt lattice saturation procedure can be geometrically 
understood in terms of a projection on the tropical linear space attached to a 
subset of the local affine Bruhat-Tits building, that we call the Gerard-Levelt 
membrane. This provides a way to compute the true Poincare rank, but also 
the Katz rank of a meromorphic connection without having to perform gauge 
transforms nor ramifications of the variable. We finally present an efficient algo- 
rithm to compute this tropical projection map, generalising Ardila's method for 
the case of the Bergman fan to the case of the tight-span of a valuated matroid. 
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Introduction 

Given a meromorphic linear differential system on the Riemann sphere, 

^- = A(z)X with A(z) e M n (C(z)) , (1) 
dz 

it is important to determine whether a singularity of A is a regular singular point for 
the system ([T]). Unlike with scalar linear differential equations, for which there is a 
purely algebraic condition on the orders of the poles of the coefficients due to L. Fuchs 
[Fuj . a system ([!} can display arbitrarily high pole orders at a regular singularity. 

Consider the differential system, expanded in the neighbourhood of the singular 
point (assumed for simplicity to be z — 0) as follows, where we put 8 = z4-, 

OX = — A i z i X with p > and A ^ if p > 0. (2) 

* E-mail : ecorel@gwdg.de 



1 



2 



The integer p is traditionally known as the Poincare rankp(A) of the system. Finding 
the type of singularity involves knowing the minimum value m(A), sometimes known 
as the true Poincare rank, of this rank under gauge transformations 

A [P] = P^AP - P^eP with P = p kZ k S GL n (K) where K = C((z)). (3) 

Several lines of research have been opened to tackle this problem. The most classi- 
cal tries to iteratively construct a suitable gauge transformation P, usually coefficient 
by coefficient in the series expansion. Featured methods rely on the linear algebra 
over C involved by equation ([3]), like Moser and continuators ( |Mol IH-Wj ). whose 
methods are widely used nowadays in computer algebra, or other researchers such as 
[Jul IBa-J-Lj . while |B-V] use Lie group theoretic tools. 

The nature of a singularity of A can also considered from point of view of mero- 
morphic connections |Dej , and especially, as a question of stability of certain lattices 
under the differential operator induced by the connection [Ml IKaj . We focus here 
specifically on the approach of saturating lattices used by Gerard and Levelt |G-L] : 
the true Poincare rank m(V) is the minimum integer k such that the sequence of 
fc-saturated lattices (recalled in section [O]) eventually stabilises. 

Recent work has shown close relations between the geometric framework of the 
Bruhat-Tits building of SL(K), for some discrete valued field such as K = C((z)), and 
tropical convexity [K-TI lJ-5-Yl IW] . We show here that Gerard and Levelt's approach 
can be formulated and efficiently computed in this framework, as a projection on the 
tropical linear space L p attached to the valuated matroid p corresponding to a given 
membrane in the Bruhat-Tits building. 

Results of Keel and Tevelev |K-T] show that, when lattices are in a same membrane 
of the Bruhat-Tits building, they are homothetic if and only if they are projected on 
the same point of a tropical linear space by an explicit nearest point projection map 
( [J-S-Yj . see also |C-G-Q[ IG-K] V We need to be able to tell whether lattices are 
actually equal. We show that the projection map respects the valuation, and projects 
only equal lattices on a same point (theorem [5]) . 

We then construct the Gerard-Levelt membrane that contains all the fc-saturated 
lattices (proposition [T]), and give a tropical version of Gerard-Levelt's lattice stabili- 
sation criterion (corollary [lj. Third, we show that the same formula amazingly allows 
for the computation of a more subtle invariant, namely the Katz rank of the connec- 
tion (theorem [3J , without having to either to compute a single gauge transformation 
nor perform the usually required ramification of the variable. 

The projection map given in | J-S-Yj has unfortunately a high complexity. We give 
in section|4]an efficient algorithm to compute this projection map onto the linear trop- 
ical space attached to a valuated matroid. We generalise the algorithmic approach 
to tropical projection developed by Ardila |Ar] for ordinary matroids, to the case of 
valuated matroids defined by Dress and collaborators |D-Wj . The algorithm presented 
here, which computes the nearest l^a -projection on the tight-span of a valuated ma- 
troid, has a wider applicability than the differential computations explained in the 
previous parts, especially in phylogenetics |D-T] . 
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1 Meromorphic connections 

A meromorphic connection is a map V : V ~ K n — > fl(V) = V ®k ^c(K) which is 
C-linear and satisfies the Leibniz rule 

V(/u) = v (g> df + fVv for / £ K and v £ V. 

The matrix Mat(V, (e)) is given by Vej = — Y%=1 e » ® f° r a basis (e). A basis 
change P £ GL„ (.fT) gauge-transforms the matrix of V by 

tt [P] = P -1 OP - P _1 dP. (4) 

Contracting with yields a differential operator V/s, and system (fTJ) is the 

expression of V_i(v) = in the basis (e). 

A lattice A in V is a free sub-O-module of rank n, that is a module of the form 

n 

A = Oei for some basis (e) of V. 

i=l 

The Poincare rank of V on the lattice A is defined as the integer 
p A (V) = -« A (A + Vo(A)) = min{fc £ N | V fe (A) c A} = max(-v(fi«) - 1, 0). (5) 

The true Poincare rank m(V) = minApA(V) characterises the nature of the sin- 
gularity of (U}, in the sense that z = is regular if and only if m(V) = 0. 

1.1 Gerard-Levelt's saturated lattices 

For any vector e £ V and any derivation r £ Der(_ftT/C), define for £ £ N the family 

Z:£(e) = (e,V T e,...,V f r e). 

The module 0^(e) induced over by Z*(e) only depends on the valuation v(t) of 
the derivation t. We can therefore restrict ourselves to the particular derivations 
Tfc = z k+1 -^ for k £ N. In this case, we put V Tfc = Vfe, and note Z|(e) and 0|(e) the 
corresponding objects. For fc ^ 1, Gerard and Levelt define the lattices 

F fe £ (A) = A + V fc A + --- + V£A. 

Note that 

n 

Pjt(A) = Zl(ei) for any basis (ei, . . . , e„) of A. 

»=i 

Theorem 1 (Gerard, Levelt). The true Poincare rank m(V) o/V is 

m(V) = min{fc e N | F£ (A) = F£ _1 (A)} for any lattice A C V. 



4 



This means that k ^ m(V) if and only if V fe (F^ _1 (A)) C Pfc _1 (A) for some 
(equivalently, any) lattice A in V, that is, that the Poincare rank on F^ _1 (A) is at 
most k. Stated otherwise, finding the true Poincare rank is finding the largest lattice 
whose Poincare rank is bounded by its index in the following sequence 

^(AjD-o^AjDA. (6) 

Let us extend this notation to multi-indices. Let £ ^ 0, and let a — (ai, . . . , at) G 
T} be an integer multi-index of length |a| = £ and weight w(a) =«! + ••• + eng. Let 
us define also the partial multi-indices au — (a\, . . . , Oj) and 

V" = V a( o--oV ai . 

Let by convention aqo = e and V e = idy for the empty sequence e. Let finally a (e) 
be the 0-module spanned by the sequence 

Z a (e) = (V Q w e) 0<J - <w . 

Lemma 1. For any a = (ai, . . . , ) € Z , one /las 

V a = z ^(«)p Q (v ) w/iere P Q (A) = X(X + cti) ■ ■ ■ (X + w(a^_i)) G Z[A]. 

Proof. The proof goes by induction on the length of the multi-index a. Let D = Vo 
The claim obviously holds for a multi- index of length 0, with P e = 1, so assume that 
there exists P Q G Z[A] such that V" = z w ^P a {D) for |a| < £ Let /3 e Z f+1 . Then 
by definition, we have 

= (w(P\ t )z w WPf) lt (D) + z w ^D o P fiit (£>)) 

= ^)(«;(^(i))+^(D)oi)) 

Indeed, P^ (D) commutes with D since it has by assumption constant coefficients. 
The result follows, since we have then Pp(X) = Pp i( (X)(X + w(f3\g)). □ 

Lemma 2. Let A be a lattice in V. For any £ G N and a G N £ , t/ie 0-module a (e) 
is spanned over fry i/ie family 

(e,z Ql V e,...,z , " (Q i f - l) V^ 1 e) . 

Proof. According to lemma[TJ the family Z^ 1 (e) is related to Z^ 1 {e) by the matrix 
P = Az Wa where 

W« = diag(0, ai,..., w(a| f _i)), 

and A is an upper triangular integer matrix with diagonal entries equal to 1, therefore 
A G SL,(Z) C GL £ (0). The families z w -Z^ x {e) and Z^e) are related by the 
matrix 

P = z~ w "P = z- w «Az w « whose entries are A tj z w{a ^- w(a ^ . 

Since A is upper triangular, and the partial sums + • • • + ctj are non-negative, the 
matrix P is in GL„(0), and therefore both families span the same O-module. □ 
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2 Tropical convexity and lattices 

Let M = {di, . . . , d m } be lines in V such that d\ + • • • + gJ to = V, and consider the 
subset of A defined by 

[M] = {t\ H 1- £ m | li is a lattice in di}. 

Following Keel and Tevelev, who call in |K-T] the set induced by [M] modulo homo- 
thety a membrane, we call this the affine membrane spanned by M. 

For a choice A — (ui, . . . , v m ) of non-zero vectors in the lines di, any lattice in 
the membrane defined by M — {di, . . . , d m } can be represented (non uniquely) by an 
integer valued point as follows: a lattice point u G Z m corresponds to the lattice 

m 

A t , = 5^0z-"i; i . (7) 

i=l 

Membranes spanned by m lines in the Bruhat-Tits building have a faithful repre- 
sentation as tropical linear spaces in m-dimensional space. 

Let (Moo =1U {oo}, ©, 0) be the tropical semialgebra, where the operations are 

x (By — min(x, y) and x y = x + y for x, y € Kqo- 

An affine membrane M and a basis (e) of V determine a valuated matroid 

p:[m] n — > Moo 

cj i — ► u(det( e ) Mo,) 

where M w = (^ 1 ,--.,u Wra ) is the subfamily of vectors of M indexed by u. To a 
valuated matroid p of rank n over [m] there is a tropical linear space L p C K.™ 
attached as follows 

L p — \ x e K™ | Vr G f ] , min p(r\{r,-}) + x n is attained twice! . (8) 

^ \n +1/ l^i^n+l J 

Depending on the authors, L p is said to be a tropical convex cone ( |C-G-Q] ) or a 
convex polytope ( |J-S-Y| ) in K™ . Both definitions mean that 

A u © /i w G L p for any A,/j £ Rqo and u, w G L p . 

According to |C-G-Q| EH EES, for id™ the formula 

Tip (a;) = min{ui G L p \ w ^ x} where the minimum is taken coordinate- wise, (9) 

defines the nearest point projection map ttl p : R™ — > L p . There are at least two 
other known ways to characterise or compute ttl p (x). 

Blue Rule. Adapting [Arj , the authors of |J-S-Y| show that hl p {%) — (wi , ■ ■ ■ , w m ) 
with 

Wi = min max(p(<7 U {i}) — p(a U {j}) + Xj). 
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Red Rule. Similarly, starting with v = (0, . . . , 0) £ R m , for every r £ Q+i) such 
that a = mini^i^„ + i p(T\{ri})+x Ti is only attained once, say at tj, compute 7 = /3— a 
where /3 is the second smallest number in that collection, and put v Ti := max(w ri ,7). 
Then ttl p (x) = x + v. 

Theorem 2 (Keel-Tevelev) . The nearest point projection map ttl p : R™ — > L p 
induces a bijection '3/ m between [M] and the lattice points in L p 

*m(A«) = vr Lp (ui, . . . ,u m ). 

In particular, M u = (z~ Ul vi, . . . , z~ Um v m ) and M u > — (z~ Ul Vi, . . . , z~ u mv m ) span 
the same lattice A if and only if ttl p (u) — ttl (u'). 

Proo f. Let A = YT=x Oz~ u 'v u and let w = ir Lp (u). According to [EE], th. 4.17 (see 
also [J-S-Y] th. 18), there exists a £ R such that Wi = v\(vi) + a for all 1 < i < m. 
By definition, v\(x) = max{fc £ Z | x £ z k A}. Accordingly, we have z~ Wi V{ £ z~ a A, 
and thus a ^ 0. By formula (J9]), we get a = and thus 

^L p (u) = (ua(Wi), ■ • ■ ,VA(v m ))- (10) 

By construction, if u' ^ u, then z~ Ui Vi = z^- Ui ~ Ui ^ z~ Ui Vi £ Oz~ Ui Vi, for 1 ^ i ^ m, 
hence A C A„/. Since in particular w ^ u holds, we get A C A w . Conversely, we 
have A w = J2"Li Qz~ Wi Vi C A. Therefore, if ir Lp (u) = ir Lp (u') then A u = A u >. The 
converse follows directly from (flOl) . □ 

3 The Gerard-Levelt membranes 

Proposition 1. Fix a basis (e) of A, and I > 0. Let [Mi] be the membrane spanned 
by the vectors (VQei)i^i^ n fi^j^e- Then Ff; (A) £ [M(\ for all k^ and I 1 si I. 

Proof. For the considered basis (e), the lattice L = Ffc (A) satisfies 

L = Q (ei) + ■ • • + Q (e„) with a = (0, k, . . . , fcl')- 

Reordering terms as (ei, . . . , e n , Voei, . . . Voe„, . . . , Vg e„), formula ([7]) and lemrnad] 
imply that L can be represented in the membrane [Mg>] by the lattice point 

(0,...,0, /c , . . . , fcj . . . , k£ . 

n times n times n times 

Since by definition, z^ VA ^v £ A holds for any v £ V, the module L can also be 
represented as an element of the membrane [Mi] by 

(0,...,0, -fc,...,-fc , . . . , -kt',. -W , u A (V£ +1 ei), va(Vo>„)). 

n times n times n times 

□ 

The lattices F^(A) for ^ I ^ n can therefore all be seen as elements of the same 
membrane [M n ]. 

Definition 1. Ma = [M n ] is called the Gerard-Levelt membrane attached to A. For 
any basis (e), the lattice F^{A) is represented by the lattice point 

u{ = (0 .,0, -fc .,-fc , . . . .,-fcl ,p A (vg +1 ei), . . .,« A (V5e„)). 

n times n times n times 
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If Mat(Vo,(e)) = A for an basis (e) of A, then Ma is described in (e) by the 
n x n(n + 1) matrix 

M=(l n A ■■■ A n ) where A k+ i = (z— + A)A k and A = I n . (11) 

The tropical projection 7ta onto the tropical linear space L\ attached to the 
Gerard-Levelt membrane Ma maps a point u to a unique representative. Check- 
ing if k ^ rn(V) requires to know if the lattice points u^ -1 and uj? represent the same 
lattice, that is 

7T A «) = TTaK- 1 ). 

Corollary 1. For any A, we have m(V) = minjfc G N | tt\(u^) = 7ta(m)! 
3.1 Tropical computation of the Katz rank 

The tropical setting is compatible with the ramification of the variable. This implies 
the following result. 

Theorem 3. Let 7ta : Ma — > L\ be the tropical nearest point projection map of the 
Gerard-Levelt membrane Ma of any lattice A onto its attached tropical linear space 
L. Then the Katz rank fc(V) of the connection V satisfies 

k(V) = min{fc 6 R + | 7Ta(u5J) = ^("fc^ 1 )} for any lattice A. 

Proof. The Katz rank is the minimum Poincare rank of the connection V# induced 
on the pure algebraic extension H — K [T)/(T N — z) of K with N = lcm(l, 2, . . . , n) 
(see e.g. [Cor] ). If we put ( for the class of T, then Mat((V h)cjl ; ( e ® 1)) = 
ATMatfV, <« , (e)). Thus if X(z) satisfies z-^X(z) = A(z)X(z) the system satisfied 
by Y(C) = d X(( N ) is 

C^Y(C)=NA(C N )Y(0. 

Put A(() = NA(( N ). The sequence (A k ) ke ® defined by relation {TTJ of iterated 
(^■-derivatives of ^4 satisfies 

A k (C)=N k A k (( N ). 

Let q be the valuatcd matroid defined by q{uj) — iy(det Af(() w ), for any n-subset u) of 
indices of the columns of M(C) with respect to the (-adic valuation w. By construction 
we have 

q(cj) = w{detM(() u ) = w(det !/((%) 
= w((detM(C) u ) N ) 
= Nw{detM{() u ) 
= Np(u). 

The lattice Njj — YaLi OHC _ "* w i ® 1 has tropical representation in L q as the pro- 
jection of the point u £ Z m with respect to the matroid g = Np. By corollary [I] 
we have m(Vfj) = min{fc G N | = 7rjv H (mJ!" 1 )}. On the other hand, 

k(V) = -^-m(Vff) holds. Therefore, we get 

k(V) = min{fc G | tt a «) = tta^- 1 )}. 

This formula holds for any extension H' of degree divisible by the denominator s of 
k(V), hence the result also holds in the limit, yielding the claimed result. □ 
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Example (Pfliigel-Barkatou). Let dX/dz = AX with 



A = 



f-5z- 2 


bz- 1 


-2Z- 1 


-9z 




3z- 2 


2z- 2 


-Az 


Az- 1 


-bz- 1 


-bz~ 2 


2 


2-2z 


-bz- 1 


Zz- 2 


-6z 



ul = (0, . . . , -3fc, -4k -4k -4k- 4k) 

v! 1 " 1 = (0,..., -3^,-6,-5,-5,-6). 

3 3 
One gets tt(u^) = tt(u^ ) k ^ — , therefore m(V) = 2 but actually k(V) = — . 



4 A projection algorithm on a tropical linear space 

The Blue and Red rules from |J-S-Y| recalled in section 13.11 have unfortunately a 
high computational complexity, since it involves loops over cardinality ("j sets. In 
our case, it is especially impractical since for the Gerard-Levelt membrane, we have 
m ~ n 2 . In this section, we present an efficient algorithm, inspired by Ardila's work 
on ordinary matroids [ Ar . to compute the projection of a point x £ R m onto the 
tropical linear space L p attached to a valuated matroid p. 



4.1 Valuated matroids 

Let us recall the setup of valuated matroids, and fix the notations that we will use. 
For the results listed in this section, we refer to |Mu-Taj . although their definition, 
following |D-Tj . comes with the opposite sign. Let E be a finite set, and a map 
p : 2 E — > Roo = RU {oo}. Let S = {B C E\p{B) / oo}. The pair (E,p) is a 
valuated matroid if 23 ^ and for B, B' S 25 and u 6 B\B' there exists v £ B'\B 
such that 

p(B) +p(B') ^ p(B U M\M) +P(B' U {u}\{v}). 

A subset B g 23 is called a basis of p. In particular, 23 is the set of bases of an 
ordinary matroid P on E, that we call the matroid underlying p. Any vector of the 
form 

X(B, v) = (p(B U - p(B), ueE) 

for some basis B and v £ E\B is a circuit of p. If A is a circuit of p, its support 

X = {e € E | X e ^ oo} 

is a circuit of the matroid P. More precisely, it is the fundamental circuit of B and v, 
that is the unique circuit of P included in B U {v}. Similarly, any vector of the form 

X*(B,v) = (p(BU{u}\{v})-p(B),ueE) 

for some basis B and v £ B is thus a cocircuit of p. 

Some important features of circuits and cocircuits of p are in fact encoded in the 
underlying matroid P. For any circuit C of P, the set of circuits of p that have C as 
support is of the form 

A + a(l,...,l) for aeR. 

Conversely, for any circuit X of p, X + a(l, . . . , 1) for a S R is a circuit of p. The 
same result applies to cocircuits. Recall the following result. 
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Lemma 3. Any circuit (resp. cocircuit) of P containing v £ E can be represented as 
the fundamental circuit (resp. cocircuit) of a basis B such that v ^ B (resp. v 6 B). 

Proof. Let C be a circuit of P. By definition, for any v s C, the set C\{v} is contained 
in some basis B. Therefore C C BU{v} holds. But there is a unique circuit satisfying 
this condition. Since the cocircuits are the circuits of the dual matroid, the same result 
holds. □ 

In what follows, we will speak by abuse of notation of the fundamental (co-)circuit 
of B and v for a valuated matroid p. This is harmless as long as the results that 
we state are invariant up to the addition of a constant. If we need to specify a 
representative, we will often use the only one with non-negative coordinates and with 
minimum coordinate equal to 0, or with some fixed value at some element of E. 

For any x <E R m , the map 

p x (B) =p(B)-J2%b 
beB 

extended to all 2 E by p x (A) — oo for A ^ 23 defines another valuated matroid on E. 

Lemma 4. If X is any circuit of p, then X + x is a circuit of p x , and if X* is a 
cocircuit of p, then X* — x is a cocircuit of p x . 

Proof. By the definition of a circuit of p, circuits of p x have coordinates 

X x {B,v) u = p x {B U {v}\{u}) - p x (B) for some Bjv 

= p(U{v}\{u}) - Xb-p(B) + J2 x b 

beBU{v}\{u} beB 

= X(B, v) u - x v + x u . 

Hence, X X (B, v) = X(B, v) + x — x v (l, . . . , 1). Similarly, we have 

X* x (B,v) u = p x (BU{u}\{v})-p x (B) 

= p(U{u}\{v})- x b -p{B) + Y J x b 

beBU{u}\{v} beB 

= X*(B, v) u + x v - x u - 

Hence, X*(B, v) = X(B, v) —x + x v (l, . . . , 1). By the projectivity property of circuits 
and cocircuits, the result is established. Since the sets of bases for p and p x coincide, 
these are indeed the only (co)circuits of p x . □ 

4.2 An algorithm for the projection on the tight-span 

A valuated matroid p : ( B ) — > Moo of rank n over a finite set E = [m] induces a 
tropical linear space L p defined by |5]) . This subspace of R™ corresponds (up to sign) 
to what Dress and Terhalle call the tight span of a valuated matroid. In this section, 
we present an efficient algorithmic method to compute the tropical projection from 
R m onto L p , that generalises results obtained by Ardila for ordinary matroids in |Arj . 

Proposition 2. Let p be a valuated matroid of rank n on [to], and let u € E. The 
following conditions are equivalent. 

i) u belongs to at least one minimal basis of p. 
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ii) u is never the unique minimum in a circuit of p. 
Hi) u is minimal in some cocircuit of p. 

Proof. = N m[) : Assume that B is a minimal basis containing u. Let C* — X*(B,u) 
be the fundamental cocircuit of B and u. By definition, we have 

C* v = p(B U {v}\{u}) - p(B) > = C*. 

That is, u is minimal in the cocircuit of B and u. 

Imp = H iiy . suppose that u in the unique minimum for p on a circuit C. Assume 
that C* is a cocircuit of p where u is minimal. By assumption, we have 

C u < C u , and C* ^ C*, for u' ^ u. 

Accordingly, C + C* has a unique minimum at u. By orthogonality of circuits and 
cocircuits f |Mu-Taj . th. 3.11, p. 204), the set of indices that minimise C + C* cannot 
have cardinality one. Therefore, the contradiction is established. 

Let us finally prove [ii]) => T7j) : consider a minimum basis B. If u £ B, let C = 
X(B,u) be the circuit generated by B and u. By assumption, the minimum in C is 
attained at v ^ u. The support of C is equal to the fundamental circuit of B and u 
for the ordinary matroid underlying p. Therefore, B U {w}\{v} is a basis of p and 

p{B U {u}\{w}) - p(B) < p(B U - p{B) for u' E C. 

Putting v = u' we get p(BU{u}\{v}) ^ p(B). Since we assumed that B was minimal, 
we get [zj. □ 

Therefore we get the following characterisation of the (finite part of the) tropical 
linear space L p . 

Proposition 3. Let x G K m , and let p be a valuated matroid of rank n on [m]. The 
following are equivalent. 

i) x G L p . 

ii) Every element of E belongs at least to one x-minimal basis of p. 

Hi) Every circuit of p contains at least 2 x-minimal elements. 

iv) Every element of E is x-minimal in at least one cocircuit of p. 

Proof. ® and lm[) are equivalent by the definition of L p (cf. |J-S-Y| ). The remaining 
assertions are obtained by applying proposition [2] to the valuated matroid p x . □ 

Note that the previous characterisation does not apply when x has an infinite 
coordinate, since p x is then no longer a valuated matroid. However, x u — oo happens 
only when u docs not belong to any basis. 

The computation of ttl v (%) can be performed independently for every element of 
the vector x. For a given u G E, there is a (unique) normalisation of a circuit C of p 
containing u such that = x u . 

Proposition 4. If u G E violates any one of the three conditions of proposition^ for 
the matroid p x , then u satisfies them for the modified vector x' — (x\, . . . , x' u , . . . , x m ) 
with 

x' = max min Cf , 

uec eec\M 

where all the circuits are normalised so that C^ — x u . Moreover, the conditions of 
proposition^ are not satisfied at u for x" — (x%, ... ,x' u — £,... , x m ) with e > 0. 
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Proof. By assumption, u is the unique x-minimum over some circuit C containing 
u. The support of such a circuit C can be defined as C — X(B,u) the fundamental 
circuit of u and a basis B $ u. The x-value at e € C of the circuit C is of the form 

= p(B U {u}\{e}) — p(B) + x e + a for some constant oel, 

so we may choose as representative of any circuit C containing u the only one such 
that C x = x u , namely the one defined by Cf = p(B U {«}\{e}) — p{B) + x e . 
Then by assumption 

x' u > min C% > CI = x u , 

e£C\{u] 

and for any circuit C B u, we have 



C? = 



Therefore 

Cf' ^ min C™ = min cf 

e£C\{u} £ eeC\{u} 

so u cannot be the unique x'-minimum of any circuit containing u. On the other 
hand, there exists a circuit C containing u such that mhi ee c\{u} C x — x' u . Putting 
x u = x 'u ~ £ f° r an y £ > 0, then u will be the x"-unique minimum over the circuit C . 
Thus x' is the smallest vector that corrects the value of i at u. □ 

Proposition 5. Let x £ K m and let B be an x-minimal basis of p. Then ui = itl p (x) 
can be computed in the following way. 




UJ; 



Xi if i G B 

Tahi u ^i{p{B U {*}\{it}) — p{B) + x u ) otherwise. 



Proof. If i € B holds, then i is x-minimal in the fundamental cocircuit X*(B,i). 
Therefore all conditions of proposition [2] apply to i. Otherwise, let X(B,i) be the 
fundamental circuit of B and i, normalised so that X(B,i)i = Xi. We have to prove 
that 

mm(X(B, i) e + x e ) — max min Cjf. 

e^i ieC e£C\{i} 

By construction, holds. Actually, it is sufficient to prove that min e5 ^(X(.B, i) e ) ^ 
min ee c'\. u\ C e for any circuit C containing i. 

Say that min^ X{B, i) e = X(B, i) u . For any circuit C containing i, there exists 
v e C\{i} such that B U is a basis. Since i B U the circuit 

C = X(B U {v}\{u},i) containing i is well defined, and C v min^^^ C e holds. 
By definition, we have 

c v = P puW\W)u{i}\H)-p(flu{i)}\M) 

= p(BU{i}\{u})-p(BU{v}\{u}) 

^ p(B U {«}\{m}) — p(B) since B is minimal 

Accordingly, we get X(B,i) u ^ C v > min eg( ^ { . } C e . □ 
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This result implies the following efficient method to compute the tropical projec- 
tion TT Lp {x) for lei" 1 . 

1. Compute a minimal basis B of p x . This can be performed by the greedy algo- 
rithm described in |D- Wj . 

2. For i 6 E\B, compute C = X(B,i). To do this, it suffices to compute the 
fundamental circuit of B and i for the underlying matroid P. 

3. Find the minimum element of C + x outside i; note that there are at most n 
non infinite elements to consider. 
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